Local density of states at polygonal boundaries of d-wave superconductors 
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Besides the well-known existence of Andreev bound states, the zero-energy local density of states 
at the boundary of a d-wave superconductor strongly depends on the boundary geometry itself. In 
this work, we examine the influence of both a simple wedge-shaped boundary geometry and a more 
complicated polygonal or faceted boundary structure on the local density of states. For a wedge- 
shaped boundary geometry, we find oscillations of the zero-energy density of states in the corner of 
the wedge, depending on the opening angle of the wedge. Furthermore, we study the influence of 
a single Abrikosov vortex situated near a boundary, which is of either macroscopic or microscopic 
roughness. 

PACS numbers: 74.45.+c,74.20.Rp,74.25.-q 



I. INTRODUCTION 

The local density of states at the boundary of a super- 
conductor is a crucial factor in many experiments, for ex- 
ample tunneling measurements. For conventional s-wave 
superconductors, the local density of states at an insu- 
lating boundary is practically the same as in the bulk. In 
particular, the specific boundary geometry is irrelevant. 
In the case of d-wave symmetry, however, the situation is 
completely different. Due to Andreev bound stateslAii, 
a drastic enhancement of the low-energy density of states 
can be observed at a straight flat surface appearing as a 
pronounced zero-bias conductance peak 4 *^. This effect 
is maximal if the d-wave nodal direction is perpendicu- 
lar to the boundary and shrinks when the orientation is 
changed 2 - 4 . For an angle of 45 degrees between nodal 
direction and boundary, the Andreev bound states dis- 
appear completely. Besides this well-known effect, it is 
important to realize that for d-wave symmetry also the 
boundary geometry itself can have strong influence on 
the local density of states. In this work we examine the 
local density of states at the surface of a d-wave super- 
conductor for some basic examples of polygonal bound- 
ary geometries and show that Andreev bound states are 
sensitive to the boundary geometry. We also consider the 
additional influence of a single Abrikosov vortex pinned 
near the boundary. 

In a previous work we have shown, that the presence 
of an Abrikosov vortex in the vicinity of the bound- 
ary of a d-wave superconductor has a drastic effect on 
the zero-energy Andreev bound states. This effect ap- 
pears as a surprisingly large shadow-like suppression of 
the zero-energy density of states between vortex and 
boundary 7 . The suppression is due to the supercur- 
rent flowing around the vortex, which locally shifts the 
quasiparticle energies and leads to a splitting of the zero- 
energy peak. A similar splitting has already been stud- 
ied for a homogeneous surface currenlp^iiSii 1 .. One of the 
purposes of the present study is to show that this vortex 
shadow effect is stable and robust even if the boundary of 
the superconductor possesses some roughness or faceting. 



The paper is divided into five sections. In the next 
section, we give a short introduction to the theoretical 
framework used for the calculations. In the third section, 
we present some zero-energy local densities of states for 
simple wedge-shaped boundary geometries and discuss 
some elementary effects. After that, we turn to the as- 
pect of boundary roughness in section four, where we con- 
sider more complicated polygonal boundary structures. 
Finally, we end with our conclusions in section five. 



II. THEORETICAL BACKGROUND 
A. The Riccati formalism 

In the limit kp£ 3> 1, the equilibrium properties of 
a superconductor are contained in the Eilenberger prop- 
agator g(f,kp,ie n ). The propagator itself can be cal- 
culated by solving a transport-like equation first intro- 
duced by Eilenberger— and independently by Larkin and 
Ovchinniko\*i^. Here, we will use the so-called Riccati 
parametrization of the Eilenberger theory, which was in- 
troduced in Ref. ^] and has proven to be very useful for 
the numerical computation of the propagator (for more 
details see Ref. 15). In order to obtain a solution for 
a given Fermi vector kp at the point r*o, the coupled 
Eilenberger equations are parametrized along a trajec- 
tory passing through ro in direction of the Fermi velocity 
corresponding to kp 



r(x) = r + x VF(kp) 



(1) 



Along this trajectory, the Eilenberger equations can be 
transformed to two decoupled differential equations of 
the Riccati type for the scalar functions a(x) and b(x): 14 

hv F d x a(x) + [2e„ + A(x)*a(x)} a(x) - A(x) = 
hv F d x b(x) - [2e„ + A(x)b{x)] b(x) + A{x)* = (2) 

Here, we have already removed the vector potential A(r) 
via a gauge transformation. This can be done for the case 
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of a single pinned Abrikosov vortex in a high-K supercon- 
ductor, since we concentrate on lengthscales smaller than 
the penetration depth A. Together with the starting val- 
ues in the bulk 



a(— oo) = 
b(+oo) = 



A(-oo) 



e„ + V4 + |A(-oo)|2 

A*(+oo) 
e n + A /4 + |A(+oo)|2 



(3) 



the functions a and b can be integrated in an easy and 
numerically stable way along the trajectory. 

The Eilenberger propagator g is related to the func- 
tions a and b by 



1 



-iri I 1 - a6 2ia 
Tab V - 2 * & -1 + ab 



(4) 



Assuming a cylindrical Fermi surface, hp depends only 
on an polar angle 9 



kp = kf 



cost/ 
sin 9 



(5) 



and the local density of states, normalized to the normal 
state density of states at the Fermi level Nq, is given by 



N(r ,E) 



f 2n d9 


~l-ab~ 


/ 




o 2^ 


1 + ab 



(6) 



Practically, this means that the integrand has to be eval- 
uated for many angles 9, each angle corresponding to one 
particular trajectory through r*o. The energy E is taken 
with respect to the Fermi level and 8 can be regarded as 
an effective scattering parameter which is proportional 
to an inverse mean free path. 

In the calculations presented below, the specific bound- 
ary geometry of the superconductor has to be taken 
into account. Since we assume that the transmission 
coefficient of all the boundaries is zero, the boundary 
conditionsiSiiLiS simplify considerably. In our case, a 
trajectory which hits the boundary is simply reflected 
like a ray of light. Thus, for the calculations in this work, 
a ray tracing procedure has been used, which allows for 
multiple specular reflections. The Riccati equations 
then have to be solved on such a multiply reflected tra- 
jectory (see for example the trajectories in Fig.[5p)). 



B. Model for the Gap function 

We consider gap functions A(r,kp) that can be sepa- 
rated into 



A(r,k F ) = A oX (kp)Tp(r) 



(7) 



Here, the symmetry function for the d x 2_ y 2-w&ve is 
x(9) = cos(2#). In the following we want to concentrate 
on the part tp(r) of the gap function, which covers the 



spatial dependence, 
and phase 



It can be factorized into modulus 



(8) 



where we call the modulus f p (r) the profile function. In 
principle, both the profile function and the phase <fi(r) 
should be calculated selfconsistently via the gap equa- 
tion. For simplicity, we have taken the profile function 
to be constant. This is a good approximation for low 
temperatures due to the Kramer-Pesch effect. Even 
for higher temperatures, the changes in the quasipar- 
ticle spectra due to this simplification are rather of a 
quantitative than a qualitative nature. Compared to the 
constant model profile function, the relevant differences 
of a more realistic profile function most often reside in 
comparatively small intervals along the trajectories. For 
the situation studied in the following, this occurs near 
the boundaries or close to the vortex center. In these 
regions, the local density of states is somewhat smeared 
out or softened, if a more realistic profile function is used. 
However, the main features to be presented below remain 
unaffected, as we have checked for specific situations. 

In the absence of a vortex, the phase <f>(r) of the 
gap function simply is zero. For the case of a pinned 
Abrikosov vortex in the system considered here, as an 
excellent approximation, the gradient of the phase is par- 
allel to the current j(r) 

W(r)||J(r) (9) 

and (f>(r) is given as a solution of the Laplace equation 

(dl + d^(r)=0 (10) 

with a "phase source" of 2tt at the vortex position. In 
other words, integrating the phase gradient along any 
closed path around the vortex position must add up to a 
total of 2ir. Additionally, <fi(f) has to fulfill von-Neumann 
boundary conditions 



d n <j>{f) = 



(11) 



since currents are not allowed to cross any of the bound- 
aries. Here, n denotes the normal vector of the boundary 
under consideration. 

To solve the Laplace equation in two dimensions, we 
use the method of conformal mapping. This analytical 
technique can be found in many mathematical textbooks, 
and in some specific physical textbooks as well. Nev- 
ertheless, in the following we give the main aspects of 
conformal mapping in short. First, the original two- 
dimensional domain, where the solution of the Laplace 
equation is sought, is regarded as a part of the complex 
plane. This domain is also referred to as physical do- 
main, and it often includes a nontrivial boundary geome- 
try with specific boundary conditions. The main idea and 
task is to map this physical domain to a simple " model" 
domain by using analytic functions. Because of the lat- 
ter, the Laplace equation itself stays invariant under this 



3 




the complex power function is analytic and that 



FIG. 1: A d a .2_ H 2-wave superconductor with wedge-shaped 
boundary geometry, a denotes the opening angle of the 
wedge. The orientation between d-wave and wedge is 
parametrized by the angle 7, which gives the rotation of the 
full-gap direction with respect to the bisecting line of the 
wedge. The origin O represents the corner of the wedge. 



mapping. Thus, if the solution of Laplace's equation with 
the correct boundary conditions is known in the model 
domain, also the original problem is solved. The solution 
in the physical domain can then be achieved by simply 
mapping each point to the model domain with the map- 
ping function and taking the corresponding value of the 
solution there. 

In our case of a wedge-shaped physical domain, the 
upper complex half space serves as the model plane and 
w, wy shall denote elements thereof. Then, the function 
$, which is given by the formula 



Wy 



w v 



w v 



w v 



is the analytical solution of Laplace's equation 



a„a„.$ = o 



(12) 



(13) 



in all points uu of the upper half plane with w ^= Wy. 
Integrating the phase gradient around the singularity at 
w = wy results in a total phase difference of 27r. Thus, 
wv can be regarded as the position of the phase vor- 
tex. Also, von-Neumann boundary conditions are satis- 
fled, since the phase gradient is parallel to the boundary 
(the real axis). Of course, formula l|12f> is analogous to 
the well-known mirror image ansatz in electrostatics, the 
first factor referring to a vortex at position wy, the sec- 
ond to a virtual antivortex at the mirrored position Wy. 

Since the correct solution of the problem in the model 
domain is given by (|12fl . the next and also last step is 
to find an analytical function which maps the physical 
domain to the upper half plane. It is well-known that 



z(w) 



(14) 



maps the upper half plane to a wedge with opening angle 
a. In particular, the positive real axis is mapped onto 
itself. Thus, for the physical domain being a wedge with 
opening angle a measured from the positive real axis, the 
sought mapping function is just the inversion 



w(z) = z~ 



(15) 



However, this mapping function is only valid for a < n. 
For opening angles it < a < 2tt one has to be careful 
about the branch cut discontinuity at the negative real 
axis. A possible mapping is then given by 



w{z) 



i(7r—a) 



(16) 



Finally, combining the solution 1)1 2|l in the model plane 
with the above mapping functions w(z), the phase factor 
of the gap function for a wedge-shaped boundary is given 
by 



tfn _ w(z) - w{z v ) w(z)* — w(z v ) 



\w(z) - w(z v )\ \w{z)* - w{z v )\ 



(17) 



On the right hand side, z = r x +ir y is the position vector 
r in complex notation. Analogous, the position of the 
vortex can be chosen with the parameter zv — tv iX + 
irv.y The opening angle a of the wedge enters via the 
mapping functions (|15|) or respectively. In principle, 
the phase factor of more than one vortex is simply the 
product of several one- vortex solutions i|17|) • However, in 
the following we will confine ourselves to the case of only 
one single vortex. 



III. WEDGE-SHAPED BOUNDARY 
GEOMETRY 

In this section, we investigate the zero-energy local 
density of states of d x 2_ a 2-wave superconductors with 
wedge-shaped boundary geometries. The situation is 
sketched in Fig. a denotes the opening angle of the 
wedge, 7 the angle between the bisecting line of the wedge 
and the direction of the maximum gap. Some examples 
for wedges with opening angles a — n/2 and a = 7r/4 
can be seen in Figs. [21 and [21 Far away from the corner, 
the zero-energy density of states at the boundaries ex- 
hibit the well-known Andreev bound states. If the nodal 
direction of the d-wave is perpendicular to the specific 
boundary line, a maximum value is reached (e.g. Fig. 
|2t), both boundaries, or[31i), upper boundary). For other 
orientations, the Andreev bound states are smaller (Fig. 
[2t>) and 03d)), or they even vanish if the direction of the 
maximum gap is perpendicular to the boundary (e.g. Fig. 

or [3^,), lower boundary). In all the calculations pre- 
sented here and in the following sections, the effective 
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FIG. 2: Zero-energy density of states of a d^.^-wave super- 
conductor for a wedge-shaped boundary with opening angle 
a — n/2. The density of states is normalized to the density 
TVo of the normal state. The effective scattering parameter 
is chosen to be 8 — O.lAo. The rotation of the d-wave with 
respect to the bisecting line is a) 7 = — tt/4, b) 7 = — 7r/8, c) 
7 = 0. 
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FIG. 3: Zero-energy density of states of a d x 2_ y 2-w&ve super- 
conductor for a wedge-shaped boundary with opening angle 
a — 7r/4. The density of states is normalized to the density 
TVo of the normal state. The effective scattering parameter 
is chosen to be 8 — O.lAo. The rotation of the d-wave with 
respect to the bisecting line is a) 7 = —n/8, b) 7 = 0, c) 

7 = 7r/8. 



scattering parameter is set to 8 = O.lAo. This leads to a 
nonzero value of the zero-energy density of states in the 
bulk and limits the absolute value of the Andreev bound 
states at the boundaries. 

It is important to realize, however, that in a range from 
the corner, which is given by the lengthscale £ = hvp/ Aq, 
the boundary geometry itself strongly influences the lo- 
cal density of states. A very counterintuitive example 
can be seen in Fig. Although the zero-energy local 
density of states reaches its maximum height along both 
boundary lines, the corner of the right-angled wedge it- 
self only exhibits the bulk value of 0.23. For the smaller 
opening angle a — tt/4, the local density of states in the 
tip of the wedge can even be lower than in the bulk: In 
Fig. Ofc), where 7 = 0, which refers to a maximum gap 
direction parallel to the bisecting line, the value is 0.11. 
This is barely larger than the bulk value corresponding 
to an s-wave superconductor. 

By just taking a look at Figs. [5] and ED one might get 
the impression, that the zero-energy density of states in 
the corner is always suppressed independently of the ro- 
tation angle 7 of the d-wave. However, this is not true 
in general. In Fig. ^we examine a wedge with opening 
angle a = 7r/3. We can clearly see that for the the nodal 
direction being parallel to the bisecting line (7 = 7r/4), 
there is a strong increase of the zero-energy local density 
of states in the corner. 



A. Local density of states in the corner of the 
wedge 

In the following, we want to concentrate on the cor- 
ner point O and systematically examine, why the local 
density of states seems to behave so strange. 

In Fig. Et), the zero-energy local density of states in 
the corner is shown. The opening angle a of the wedge- 
shaped boundary geometry corresponds to the vertical 
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FIG. 4: Zero-energy density of states of a d x 2_ y 2-w&ve super- 
conductor for a wedge-shaped boundary with opening angle 
a — 7r/3. The density of states is normalized to the density 
TVo of the normal state. The effective scattering parameter is 
chosen to be 8 = O.lAo. The rotation of the d-wave with re- 
spect to the bisecting line is a) 7 — 0, b) 7 = n/8, c) 7 — 7T/4. 



axis. The rotation angle 7 of the d-wave is varied from 
to 7r/2 along the horizontal axis. The absolute max- 
imum value is right in the center of the figure. Here, 
a = 7r and 7 = 7r/4 correspond to the well-known maxi- 
mum Andreev bound states at a straight boundary line. 
The most interesting part of the picture can be seen at 
smaller opening angles a < n. If the nodal direction 
of the d-wave is kept fixed parallel to the bisecting line 
of the wedge (i.e. 7 = 7r/4) and the opening angle a 
is reduced, then the zero-energy local density of states 
clearly oscillates between high maxima and minima. A 
cut along this line showing the oscillations in the range 
0<a<7r/2is shown in Fig. Clearly, the minima of 
the oscillating zero-energy local density of states in the 
corner appear for opening angles a+ = 7r/(2n), while the 
maxima are very close to the angles a~ = 7r/(2n — 1). 
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FIG. 5: a) Zero-energy density of states in the corner O of 
a d x 2_ y 2-wave superconductor with wedge-shaped boundary 
geometry, a denotes the opening angle of the wedge, 7 the 
rotation of the maximum gap direction with respect to the 
bisecting line of the wedge. The global maximum in the mid- 
dle of the picture corresponds to the well-known maximum of 
the Andreev bound states at a straight boundary line, with 
the nodal direction of the d-wave being perpendicular to the 
boundary, b) Two examples for the relations I18H and 119> . 
The upper panel shows, that for an opening angle of a — n/2 
an incoming trajectory always leaves the corner at the same 
angle. For an opening angle of a — 7r/3, the directions of 
the incoming and outgoing trajectories are symmetric with 
respect to the bisecting line of the wedge (lower panel). 



B. Explanation of the oscillations 

For opening angles a < it, multiple reflections of a 
quasiparticle trajectory occur, because the boundaries 
are hit several times. Generally, for a given opening an- 
gle a, the exact path of a quasiparticle trajectory has to 
be determined by a raytracing algorithm. However, there 
are some specific opening angles a n — n/n, which have 
useful properties and are the key in understanding the 
oscillations. We denote 4> S [0: a ] the angle of an incom- 
ing trajectory before any reflection has occured, and call 
tp G [0, a] the angle of the final outgoing trajectory after 
all the multiple reflections. Then, the specific opening 
angles a+ = 7r/(2n) have the property, that always 



(18) 



while for a n = 7r/(2n — 1) we always have 

%j} = a-(j) (19) 

In the first case, incoming and outgoing trajectories are 
parallel. This means that the gap value is the same on the 
incoming and outgoing part of the trajectory. In contrast, 
in the second case the angles of incoming and outgoing 
trajectories are symmetric with respect to the bisecting 
line of the wedge. For 7 = 7r/4 (nodal direction along the 
bisecting line) this means, that the gap possesses opposite 




FIG. 6: Oscillation of the zero-energy density of states in 
the corner O of a d a .2_ H 2-wave superconductor with wedge- 
shaped boundary geometry. The nodal direction of the d- 
wave is fixed parallel to the bisecting line of the wedge, i.e. 
7 = 7r/4. Decreasing the opening angle a leads to oscillations 
of the density of states in the corner, with minima appearing 
at at = 7v/(2n) and maxima close to = 7r/(2n — 1). Both 
upper and lower envelope (dotted lines) are given analytically. 



sign on the incoming and outgoing part of the trajectory, 
which results in large contributions to the zero-energy 
density of states in the corner. Two typical examples 
for the different symmetry behaviour of in- and outgoing 
trajectory directions are shown in Fig. |5J>)- 

Now we turn to the calculation of the local density 
of states in the corner point O of the wedge. For the 
trajectories passing through this point, it is an excellent 
assumption to neglect all the complicated details of the 
multiple reflections in the corner completely, since they 
occur on a very small lengthscale along the trajectory 
path. The only important things to keep are the rela- 
tions between the angles of incoming and finally outgoing 
trajectories. Thus, the effect of the wedge-shaped bound- 
ary geometry is only a mixing of the bulk values © of 
a and b belonging to different angles. If we concentrate 
on the specific opening angles a n with the relations H18|) 
and i|19|) . it is possible to derive an analytic expression 
for the local density of states in the corner. 

The oscillations in Fig. H3 appear for the d-wave orien- 
tation 7 = 7r/4, where the nodal direction is parallel to 
the bisecting line of the wedge. Based on relation (|18H for 
opening angles a+ = 7r/(2n), we find the lower envelope 



1 r /2 . . 1 

a 



(20) 



-a/2 



A 



-F(a,-A%/6 2 



1 + sin^ 2<j> 



Analogous, the upper envelope, which connects the val- 
ues of the local density of states for opening angles 



a n = 7r/(2n — 1), is given by 



^Corne r (£ = 0,7 = VT/4) 



(21) 



/2 

i 

-a/2 



A 2 

1 + -=£ sin 2 2cj> 



-£(a,-A 2 /^ 



Here, F and -E are Elliptic Integrals of the First and Sec- 
ond Kind. Both upper and lower envelope are shown as 
the dotted lines in Fig. HJ1 As already mentioned before, 
we chose the effective scattering parameter S — 0.1A . 
For the lower envelope, all the bulk values stemming from 
an interval of length a around the nodal line contribute 
to the local density of states in the corner. Because of 
that, for a = tt/2 the whole quasiparticle spectrum in 
the corner is the same as in the bulk, since contribu- 
tions from the full cJ-wave are included. In particular, 
the zero-energy density of states in the corner is just the 
bulk value. When the opening angle a is reduced, how- 
ever, the result approaches the normal state value 1. The 
upper envelope is dominated by bound states, which oc- 
cur by the mixing due to relation l|19f) . since the gap 
function changes sign at the bisecting line of the wedge. 
The highest contribution to the bound states is confined 
by the opening angle a. Thus, after the maximum for 
a = 7r at the flat boundary, the upper envelope shrinks 
to the normal state value 1 for smaller opening angles. 

If the maximum gap direction of the ci-wave is parallel 
to the bisecting line of the wedge, i.e. 7 = 0, the sign of 
the gap function is symmetric with respect to the bisect- 
ing line. Then, the mixing i|19fl of the specific opening 
angles a~ generates nothing but bulk contributions, too. 
Consequently, for all the opening angles a n — n/n, no 
bound states appear at all. The corresponding values of 
the zero-energy local densities of states in the corner are 
on the same curve for both at and olZ: 



iV C omer(£ = 0,7 = 0) = - [ ' 

a J - a 



a/2 

t 

/2 



1 



1 + =$■ cos 2 2<j> 
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An 



-F(a,l/(1 + <5VA^)) 
ii , i a 



(22) 



For ol\ — 7r and again for a-i — tt/2, we have the 
same density of states both in the bulk and in the cor- 
ner. For smaller opening angles, only the contributions 
from the full-gap direction remain. The local density of 
states in the corner approaches the shape of a bulk s- 
wave spectrum. Because of that, the zero-energy weight 
approaches the s-wave bulk value, as already pointed out 
before. Although there are no zero-energy bound states 
in the corner of the wedge for the specific opening an- 
gles a n = tt/h, multiple reflections lead to considerable 
bound states for angles tt/2 < a < it. For arbitrary 
opening angles a < tt/2, however, the deviation of the 
zero-energy density of states in the corner from Eq. 1L'2I) 
determined by the a n is small. 

To summarize this section: Whether Andreev bound 
states appear in the corner of a wedge or not depends 




FIG. 7: Zero-energy density of states of d x 2_ y 2-wa.ve super- 
conductors with wedge-shaped boundary geometries. A sin- 
gle Abrikosov vortex is pinned at a distance of 2 £ from each 
boundary line. All pictures presented here correspond to a 
specific wedge and orientation already shown earlier without 
a vortex. The same color palettes have been used for both the 
corresponding plot without vortex and the plot given here, a) 
This picture refers to Fig. [2J:) . The vortex casts a shadow-like 
suppression on the Andreev bound states at each boundary, 
b) This picture refers to Fig. 0Ji). Again, the Andreev bound 
states at the boundaries are locally suppressed, c) This pic- 
ture refers to Fig. The bound states in the corner in- 
duced by the boundary geometry are also affected but remain 
still considerable. Please note, that b) and c) are given in a 
smaller scale than the corresponding pictures. 



on both opening angle a of the wedge and orientation 
7 of the d-wave. For a < tt/2, Andreev bound states 
in the corner are suppressed in most of the parameter 
space. However, near the specific opening angles a~ = 
tt /(2n — 1) and for orientations about 7 = 7r/4, Andreev 
bound states get induced in the corner. Fig.[S^) provides 
a map, showing for which combinations of angles Andreev 
bound states appear. 



C. Wedge and Vortex 

We already discussed the influence of a single 
Abrikosov vortex, which is pinned near a straight smooth 
boundary, on the local quasiparticle spectrum in a previ- 
ous worki. The result is a suppression of the zero-energy 
Andreev bound states in a shadow-like region extending 
from the vortex to the boundary. This effect is due to the 
flow field of the phase gradient around the vortex, which 
leads to a local shift of the quasiparticle energy along 
the trajectory. As a consequence, the zero-energy An- 
dreev bound states at the boundary are suppressed and 
the spectral weight is shifted towards higher quasiparticle 
energies as discussed in Ref. 0- 

Some examples of a single Abrikosov vortex, which 
is pinned near a wedge-shaped boundary geometry, are 
shown in Fig. [7| Here, we have calculated the phase 
distribution around the vortex using the conformal map- 
ping procedure described in section II B. Each picture 
corresponds to a picture with the same wedge geometry 
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and d-wave orientation already presented in earlier fig- 
ures. As a main effect, the existence of the vortex and 
its inhomogeneous phase lead to a local suppression of 
the zero-energy Andreev bound states at the boundaries. 
The center of the suppression seems to be approximately 
that point of the boundary, which lies closest to the vor- 
tex. 

If we concentrate on the very corner of the wedge, how- 
ever, the characteristic changes of the local density of 
states induced by the boundary geometry itself are only 
slightly affected by the presence of the vortex. Of course, 
the induced bound states in the corner of Fig. for 
example, get reduced by about 30 % in Fig. 0;) because 
of the vortex. Nevertheless, the zero-energy density of 
states is still strongly increased in the corner compared 
to the bound states at the boundary nearby. An induced 
suppression of the zero-energy density of states in the 
corner is practically not affected by the vortex as can 
be seen in Figs. Ek) and b). It is important to realize, 
however, that the shown suppressions of the zero-energy 
local density of states are not of the same nature. On the 
one hand, the suppression due to the vortex is the result 
of a splitting of the sharp zero-energy peak in the quasi- 
particle spectrum. On the other hand, the low value of 
the zero-energy density of states because of the boundary 
geometry is just a natural consequence, when the quasi- 
particle spectrum in the corner is given by a d- or even 
s-wave bulk spectrum. 

Some further examples of a wedge with opening angle 
a = 3ir/2 and an Abrikosov vortex at different positions 
can be seen in Fig. [8] In the top row, the shadow effect 
can be observed again. Below, there are no bound states 
at the boundaries at all, and the influence of the vortex 
rather leads to a very small increase of the zero-energy 
density of states. 



IV. POLYGONAL BOUNDARY STRUCTURE 

In this section, we want to study the influence of sur- 
face roughness and surface faceting on the vortex shadow 
effect introduced in Ref . in order to see how stable this 
effect is. There have been several different suggestions in 
the literature how surface roughness can be implemented 
within Eilcnberger theory (see for example Refs. 8 18.20). 
Here, we focus on the zero-energy density of states in the 
vicinity of a boundary line for two different models of ei- 
ther microscopic or macroscopic roughness with respect 
to the lengthscale of the coherence length. We begin with 
a very simple model of microscopic roughness in the next 
subsection. After that, we present results for a polygonal 
boundary structure, which can be regarded as a macro- 
scopically rough, faceted surface. In contrast to previ- 
ous work we take into account the complicated current 
redistribution along the faceted surface using conformal 
mapping techniques. 



a) b) c) 




FIG. 8: Zero-energy density of states of a d x 2_ y 2-w&ve super- 
conductor for a wedge-shaped boundary with opening angle 
a — 37r/2. A single Abrikosov vortex is situated at a vertical 
distance of 2£ from the corner of the wedge. The horizontal 
distance of the vortex is a) 2£, b) £, c) 0. In the top row, the 
nodal direction of the d-wave is perpendicular to both bound- 
aries, which corresponds to 7 = 0. Below, the nodal direction 
is parallel to the bisecting line, i.e. 7 = n/4. 



A. Microscopic roughness 

A simple model of a microscopically rough surface is 
given by an arrangement of randomly tilted tiny mirrors 
along the boundary line. The size of each mirror is taken 
to be about £/40, which is of an atomic lengthscale for 
high-T c superconductors. Each mirror is randomly tilted. 
We have chosen a Gaussian distribution between — 7r/4 
and +7r/4 for all the tilt angles, with a tilt angle of cor- 
responding to a parallel alignment of mirror and bound- 
ary line. Additionally, we made the simplification, that 
the mirrors do not extent into the superconductor. Their 
only function is to reflect an incoming trajectory accord- 
ing to the mirror orientation, thus producing a kind of 
random scattering at the boundary. There are some fur- 
ther consequences of the simplification made above. On 
the one hand, the flow field inside the superconductor 
is not affected by the boundary roughness and is sim- 
ply the one of a wedge-shaped boundary geometry with 
opening angle a — ir. On the other hand, multiple re- 
flections at the boundary and the corresponding effects 
presented in the previous section are excluded. However, 
it is important to note, that for any boundary orienta- 
tion, there exist zero-energy bound states in some regions 
of the boundary line because of the microscopic rough- 
ness. This is in agreement with earlier calculations of the 
zero-bias conductance peak for a rough surface^. 

How is the shadow effect due to a vortex affected by 
such a microscopically rough surface? The zero-energy 
local density of states for a boundary line with one spe- 
cific randomly generated mirror arrangement and an ad- 
jacent Abrikosov vortex is shown in Fig. |5J In col- 
umn a), the nodal line of the d-wave is perpendicular 
to the boundary, which corresponds to a (110) surface or 
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FIG. 9: Zero-energy density of states of a d x 2 _ y 2-wave super- 
conductor in the vicinity of a surface. This corresponds to 
an opening angle of the wedge of a — n. A single Abrikosov 
vortex is located at a distance of 2£ from the boundary. In 
the top row, the boundary is smooth. In the lower row, the 
boundary is microscopically rough, a) The nodal direction 
of the d-wave is perpendicular to the boundary, which cor- 
responds to 7 = 7r/4. b) The rotation angle is 7 = tt/6. 
Compared to a), the orientation changed by 15 degrees. 




conformal mapping 



FIG. 10: The conformal mapping procedure used to ob- 
tain the phase distribution around an Abrikosov vortex in 
the vicinity of a boundary. 
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7 = 7r/4, respectively. Although the appearance of the 
Andreev bound states along the boundary line changes 
in a characteristic way because of the microscopic rough- 
ness, the suppression of the Andreev bound states in a 
shadow- like region clearly persists. In b), the orientation 
between boundary and d-wave is changed by 15 degrees. 
The suppression of the bound states remains, but it is ob- 
viously no more symmetric, neither for the smooth nor 
the microscopically rough boundary. Since the modulus 
of the phase gradient around the vortex is symmetric, 
this asymmetry cannot be explained by a Doppler-shifted 
energy spectrum, where only local surface currents are 
taken into account. In fact, generally the whole quasi- 
particle "history" along its trajectory is very important 
and may not be neglected. Only in the case of a homoge- 
neous flow fieloS, for example far away from the vortex, 



FIG. 11: Zero-energy density of states of a d x 2_ y 2-wave su- 
perconductor with a polygonal boundary geometry, which can 
be regarded as a boundary line with macroscopic roughness 
(faceting), a) No vortex is present, b) A single Abrikosov 
vortex is situated 2£ away from the averaged boundary line. 



it becomes sufficient to consider solely the local current. 

B. Macroscopic roughness 

We now examine the zero-energy local density of 
states of a (i-wave superconductor along a macroscopi- 
cally rough boundary line. As an adequate model, we 
take a superconductor with a polygonal boundary ge- 
ometry, consisting of piecewise smooth facets, each of a 
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length of the order of £. In contrast to the model of 
a microscopically rough surface presented above, the ex- 
tension of the polygonal boundary geometry is taken into 
account completely. For the more simple wedge-shaped 
boundary in section III, we obtained the phase of the 
gap function in an analytical way by conformal map- 
ping (cf. section II B). For the more general case of a 
polygonal boundary, it becomes necessary to calculate 
the Schwarz-Christoffel mapping 21 numerically. We used 
an implementation by Driscoll 2 ^ to obtain the correct 
corresponding phase. This toolbox uses the analytically 
known solution of the phase distribution for a vortex sit- 
ting in the center of a disk and conformally maps this 
solution onto the vortex close to the polygonal bound- 
ary, as illustrated in Fig. ED For the disk, the lines of 
constant phase are just radial lines, and the phase gradi- 
ent corresponds to concentric circles, fulfilling the bound- 
ary condition that the phase gradient is parallel to the 
surface. The toolbox provides the mapping which maps 
the center of the disk onto the vortex position and the 
boundary of the disk onto the polygonal boundary. The 
resulting lines of constant phase and phase gradient af- 
ter this mapping are shown in Fig. I1UI right panel. This 
means, that within the approximation @, the currents 
due to a single Abrikosov vortex are forced to be parallel 
to each facet of the polygonal boundary. Furthermore, 
the polygonal boundary allows for multiple reflections of 
quasiparticle trajectories. Thus, the effects discussed in 
section III can be found here as well. 

In Fig. Illfa l and b), the zero-energy local density 
of states of the d-wave superconductor with polygonal 
boundary structure is shown. The main direction of the 
boundary is oriented in such a way, that facets parallel 
to it exhibit maximum Andreev bound states. If we com- 
pare the local density of states shown in Fig. Illb). where 
no vortex is present, to the case with Abrikosov vortex 
presented in b), many of the effects already discussed 
above can be found again (cf. for example Fig. 0). The 
energy spectrum in the corner of a right-angled wedge is, 
independent of its orientation, the d-wave bulk spectrum. 
Thus, the nearly right-angled parts of the boundary ex- 
hibit a very low zero-energy density of states, which is 
hardly affected by the presence of the vortex. The main 
effect of the vortex is again a strong suppression of zero- 
energy bound states in its shadow region. This is because 
the zero-energy spectral weight is effectively shifted to- 
wards higher energies. Bound states at the boundary, 
which are further away, are reduced more moderately. 

V. CONCLUSIONS 

In this work, we investigated the influence of both the 
boundary geometry and a single pinned Abrikosov vor- 
tex on the zero-energy local density of states along the 
boundary of a d^^^-wave superconductor. We found 
that a wedge-shaped boundary can induce a quasiparti- 
cle spectrum in the corner of the wedge, which is com- 



pletely different to that farther away from the corner. 
The main mechanism for this effect is the multiple reflec- 
tion of quasiparticle trajectories. In the corner, this can 
effectively be described as a mixing of bulk trajectories 
belonging to different angles. Depending on the opening 
angle of the wedge and the orientation between d-wave 
and boundary geometry, this mixing can either lead to 
bulk spectra in the corner (and thus to an absence of 
bound states) or to the presence of induced zero-energy 
bound states. In Fig.[S]i), we have provided a map, show- 
ing for which combinations of opening angle and d-wave 
orientation bound states appear. If the wedge is oriented 
in such a way, that the nodal direction of the d-wave 
is parallel to the bisecting line, the zero-energy density 
of states in the corner strongly oscillates as a function 
of the opening angle of the wedge. Although all results 
presented here have been obtained using a model gap 
function, a fully selfconsistent calculation is not going to 
change the results qualitatively. The amplitude of the os- 
cillations mentioned above, for example, should be some- 
what reduced, especially for smaller opening angles. But 
in any case, it should be possible to observe the oscilla- 
tions at least for larger opening angles. 

Another example of an induced quasiparticle spectrum 
is a wedge with opening angle n/2: The local density of 
states in the corner of a right-angled wedge is always 
given by a d-wave bulk spectrum, even if the orientation 
allows the highest Andreev bound states at both bound- 
ary lines. In addition to this geometrically induced pro- 
hibition of bound states, we presented the influence of 
a single pinned Abrikosov vortex on the zero-energy lo- 
cal density of states at the boundary. Bound states at 
the boundary are locally suppressed, in particular in a 
shadow-like region close to the vortex. However, this 
kind of suppression is due to a splitting of the high zero- 
energy peak of the bound states. Thus it is of a different 
physical origin. The spectral weight is effectively shifted 
towards higher energies, because the quasiparticles "see" 
a locally varying energy shift along their trajectory, which 
stems from the locally varying flow field around the vor- 
tex. Both effects can be observed independently in our 
results. 

We have studied the influence of two types of surface 
roughness on the vortex shadow effect: microscopic, ran- 
dom scattering at the surface as well as a more macro- 
scopic faceting of the boundary. In the second case we 
have taken care to calculate the flow field with the appro- 
priate boundary conditions using the conformal mapping 
technique. In both cases the vortex shadow effect clearly 
persists. 
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